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Abstract 

The equations of the restricted three-body problem describe the motion of a massless particle 
under the influence of two primaries of masses 1 — fi and ix, < /i < 1/2, that circle each other 
with period equal to 2ir. When /z = 0, the problem admits orbits for the massless particle that 
are ellipses of eccentricity e with the primary of mass 1 located at one of the focii. If the period 
is a rational multiple of 2-7T, denoted 2np/q, some of these orbits perturb to periodic motions for 
fi > 0. For typical values of e and p/q, two resonant periodic motions are obtained for fi > 0. 
We show that the characteristic multipliers of both these motions are given by expressions of 
the form 1 ± ^/C(e,p, q)p, + 0(/i) in the limit fi — > 0. The coefficient C(e,p, q) is analytic in e 
at e = and C(e,p,q) = 0(e' p ~ q \). The coefficients in front of e' p_g ', obtained when C(e,p,q) 
is expanded in powers of e for the two resonant periodic motions, sum to zero. Typically, if one 
of the two resonant periodic motions is of elliptic type the other is of hyperbolic type. We give 
similar results for retrograde periodic motions and discuss periodic motions that nearly collide 
with the primary of mass 1 — fx. 

1 Introduction 

Resonance between the periods of planets, satellites, and other objects is a feature of celestial 
mechanics. Among the many extra-solar planets discovered recently, many types of resonance have 
been observed. Here we consider resonant periodic motions in the restricted three-body problem, 
with the principal aim of gaining an understanding of the instabilities caused by resonance. 

The restricted three-body problem is about the motion of a massless particle under the influence 
of primaries of masses 1 — fi and \i that circle each other at a distance equal to 1 and with period 2tt. 
We will assume < jx < 1/2 and in much of the discussion /x is small and positive. Instead of the 
inertial frame of reference, it is often convenient to use a frame of reference which rotates with the 
two primaries and which is centered at the center of mass of the two primaries. The Hamiltonian 
for the motion of the massless particle in a rotating Cartesian frame is 

H = ^(Pl + P 2 y) + VPx - XPy ~ ~ (1.1) 



*Department of Mathematics, University of Michigan, 530 Church Street, Ann Arbor, MI 48109, U.S.A., 
divakar@umich.edii. This work was supported by the N.S.F. and by a research fellowship from the Sloan Foun- 
dation. 



1 



where Ao = ((x + [i) 2 + y 2 ) 1 ^ 2 and A\ = ((x — 1 + fi) 2 + y 2 ) 1 / 2 . The generalized momenta 
corresponding to x and y are p x and p y . The primaries of mass 1 — fj, and ^ are located at (— /i, 0) 
and (1 — ^,0), respectively, and Ao and Ai are the distances from the massless particle to the 
two primaries. We consider the motion only in regions where the distances Ao and Ai are both 
bounded away from zero. Therefore (1 — //)/Ao + /VAi can be expanded in powers of fj, as 

+ fj,n'(x, y, fi) = - = ■ —j, + fjSl(x, y) + 0(/u 2 ) (1.2) 



{x 2 + y 2 y/ 2 r V ""^ ( x 2 +y 2)l/2 
where 



((X- l) 2 +y 2 )V2 ( x 2 + y 2)3/2 ( x 2 +y 2)l/2- 

The change of variables p x <— — p x , x <— x, p y <— p y , y <— — y is symplectic with multiplier — 1 
and leaves the Hamiltonian H unchanged. Therefore if x(t) = X(t),y(t) = Y(t), < t < T, is a 
solution of Hamilton's equations formed using H, then x(t) = X(—t),y(t) = —Y(—t), —T <t<0, 
is also a solution. 

If r and 9 are polar coordinates then x = rcos9, y = rsin#, and the Hamiltonian becomes 
H = \ ( R2 + ^ ) " ° ~ \ ~ ^ (r ' 6) + ° ( ^ 2) ' 

where 

^, 1 cos^ 1 1 cos 6 1 

) = -r 2 = 77~T~2 n mT/2 - — • ( L3 ) 

Ai r z r (1 + r z — 2r cos y) i ' z r z r 

The discrete symmetry of this Hamiltonian is given by R <— — i?,, r <— r, G <— G, and <— — 0. 
This transformation too is symplectic with multiplier —1. 

Given R,G,r,9, another set of variables better suited for perturbative calculations can be 
defined as follows. Let Hq = (R 2 + G 2 /r 2 )/2 — 1/r. Define the variables a,e,E,u,g using the 
equations 

e = y/l + 2G 2 H 
o(l-e) =G 2 /(l+e) 

r = a(l — e cos E) 
cos E — e 



cos z/ 



1 — e cos -E 

9 = 0- v. (1-4) 

When the third and fourth equations in H1.4j) are solved for i£ and z/, there are typically two 
solutions in the interval [0, 2tt). Exactly one of these can be chosen by heeding the signs of G 
and R. Until the last section, we assume < e < 1. The variables a, e,u, g can be used instead 
of R,G,r,6. When /i = 0, the orbit (of the massless particle) is an ellipse in the inertial frame 
with eccentricity e and semimajor axis a; the angles E and v are the eccentric and true anomalies, 
respectively; and g is the argument of the perihelion in the rotating frame. 

The change of variables from R, G, r, 6 to a, e, u, g, in the region of phase space with < e < 1, 
is not symplectic. To obtain a symplectic change of variables in that region, define the variables 
L, I using the equations 

L = ±v^ I = E - esmE. (1.5) 
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Figure 1: The dotted line above is a circle of radius 1. 

As further explained in Section 3, the sign of L is chosen to be the same as that of G, the angular 
momentum. The change to L,G,l,g is symplectic and the Hamiltonian becomes 



H 



2L 2 



G - fiQ(L, G, I, g) + 0(p 2 



(1.6) 



Simple expressions for the perturbation Q in terms of x, y and in terms of r, 9 have been given in 
(|1.2|) and IJl.Mjl . In Sections 2 and 3, we partially develop f2 as a trigonometric series in I and g. 
The variable I is called the mean anomaly. The discrete symmetry of the Hamiltonian Q1.6JI is given 
by L <— L, I < Z, G <— G, g < g. This transformation too is symplectic with multiplier —1. 

Let n = 0. If L = (p/q) 1 / 3 , where p and q are relatively prime positive integers, or equivalently 
a = (p/q) 2 / 3 , and 0<e<lorG = (a(l — e 2 )) 1 / 2 , then the massless particle moves on an ellipse in 
the counterclockwise sense. In the rotating frame, this ellipse itself rotates in the clockwise sense 
with angular velocity equal to 1 as shown in Figure ^ The orbit has a period equal to 2np in the 
rotating frame as the period of motion on the ellipse is in resonance with the period of rotation 
of the frame. Arenstorf and Barrar [3] proved that the periodic orbit persists for fj, > and [i 
small enough, if I and g are multiples of tt at t = 0, and if the orbit for fx = does not pass through 
the location of the primary of mass //. 

Since I and g are angular variables, there are four possibilities for I and g at t = 0, namely 
I = ra;7r and g = n g ir with n\ and n g either or 1. Only two of these are distinct. They are depicted 
on the Poincare section g = in Figure [2j If p is odd the two distinct orbits are obtained with 
n g = and n\ either or 1. When p is even, we can take n\ = and n g either or 1. Over its 
period, the orbit intersects the Poincare section exactly p times. The plots in Figure correspond 
to p = 1 and p = 2. 

We investigate the linear stability of resonant periodic motions, such as the ones depicted in 
Figure [21 for fi > 0. The existence and stability of the nearly circular orbits, which are related to 
the nearly elliptic ones considered here, is discussed in [3], |Hj, [0], and [TT] . 

In Section 2, we prove that the characteristic multipliers of the resonant periodic orbits are 



3 




Figure 2: Periodic points on the L-l plane with p = 1 and p = 2, respectively. 

given by the expression 1 ± ^C(e,p, q)fi + O(fi) in the limit fi — > where 

6vrp 5 / 3 f 2n Pd 2 n(r,9) , 

In the integral above, expression Q1.3[) for f2(r, 9) must be used to evaluate the partial derivatives 
of S7(r, 0) with respect to 6. After that, the variables r and 6 must be expressed in terms of 
L,G,l,g using Q1.4JI and Q1.5JI . The dependence of I and g on i is obtained from (|1.6|) with /j, = 

as /(t) = gi/p + n/7r and g(t) = — t + n g ir, and L = {p/q) 1 ^ 3 and G = {p/q) 1 ^ i {l - e 2 ) 1 / 2 do not 
vary with t. 

If p/q 1/1, C(e,p,q) is analytic at e = 0. In Section 2, we prove that all powers of e with 
exponent less than \p — q\ have coefficients equal to in the power series expansion of C(e,p,q) 
about e = 0. The coefficient of e' p ~ 9 ' can be obtained as a finite sum of hypergeometric terms if 
q > 1. The coefficients of e' p-<? ! for the two distinct types determined by the choice of ni and n g 
sum to zero. It seems to be the typical case that if one of the two types of periodic motions is 
elliptic then the other is hyperbolic for small p. This is depicted in Figure El by the use of crosses 
for hyperbolic points and circles for elliptic points. In later work, we will demonstrate the existence 
of homoclinic points in the resonance band around L = (p/q) 1 / 3 - 

In Section 3, we give similar results for retrograde periodic motions. For retrograde periodic 
motions, C(e,p,q) = 0(e p+q ) as e — > 0. In Section 4, we consider periodic motions near collision 
with the primary of mass 1 — [i. The variables L,G,l,g are not a valid choice of variables at 
collision. Giacaglia |Hj solved the Hamilton-Jacobi equation for the regularized Hamiltonian with 
[i = and derived action-angle variables that are valid near collision. The existence of resonant 
periodic motions near collision was proved by Schmidt ^0]- We solve the Hamilton-Jacobi equation 
again, using the geometric approach of Arnold 2,, to correct a calculation and to put it on a more 
secure basis. We discuss the existence of periodic motions that collide with the primary of mass 
I- p. 
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2 Direct periodic motions 



Let n = 0. The Hamilton's equations formed using Ql.fij) are / = L~ 3 , g = — 1, L = 0, and G = 0. 
Consider the solution with the initial conditions L(0) = (p/q) 1 / 3 , G(0) = (p/q) 1 / 3 ^ — e 2 ) 1 / 2 , 
^(0) = run, and g(0) = n 5 7r, where p and (? are relatively prime positive integers, < e < 1, and 
ni,n g are either or 1. This solution has a period equal to 2np. Assume that Ai > along this 
orbit. Then O' in (|1.2|) is uniformly bounded in a neighborhood of the orbit for ii small enough. 
Therefore, for initial conditions chosen in a neighborhood of this solution and for a finite interval 
of time, the solution is analytic in [i at fx = 0. The assumption Ai > is valid for p/q ^ 1/1 if e 
is small enough. For given p/q, ri[, and n g , it is valid for all except finitely many values of e PQ. 

The solution is unchanged by the discrete symmetry L <— L, I < I, G <— G, g < g, at both 

t = and t = up. For ii > and [i sufficiently small, it is possible to perturb L(0) and the time of 
flight irp such that the solution is unchanged by the discrete symmetry at t = and t = irp + O(fi) 
[2]. The existence of a unique family of periodic solutions for fi > and /i sufficiently small with 
period T = 2itp + 0(p) and initial conditions L(0) = (p/q) 1/3 + 0{jt), G(0) = (p/q) 1/3 (l - e 2 ) 1 / 2 , 
1(0) = niTT, and g(0) = n g ir follows. 

2.1 Characteristic multipliers 

To obtain the characteristic multipliers of this family of periodic solutions, we consider the return 
map to the Poincare section g = n g Tr. The Hamilton's equations that correspond to (|1.6j) are 

i = l-' 3 - fin L + o(fi 2 ), g = -i-fin G + o(^), t = ^ l + o{^ 2 ), g = ^o, + o(^ 2 ). (2.1) 

Since H is a first integral of these equations, we can fix H and solve (|1,6|) for G using the implicit 
function theorem and identify the Poincare section with the L-l plane. Figure depicts the inter- 
sections of periodic solutions with p = 1 and p = 2 with the Poincare section. We refer to these 
intersections as periodic points. Since g = — 1 for [i = 0, the implicit function theorem implies the 
existence of the return map over large regions of the L-l plane and certainly in a neighborhood of 
the periodic points. The Poincare map leaves the area element dl dL invariant. When [i = 0, the 
return time is 2ir for all points in the L-l plane. 

When /i = 0, the pth return map maps the periodic point on the L-l section with L(0) = (p/q) 1 / 3 
and 1(0) = niir to itself. The pth return time denoted by T is 2np. By (f2~Tj) . l(t) = 1(0) + tL(0)~ 3 
and L(t) = L(0). Therefore, 

d(L(T),l(T)) ( 1 0\ 

d(L(0 )m = [-*ffi lj M 

at L(0) = (p/q) 1 / 3 , 1(0) = niir, which implies that the characteristic multipliers of the pth return 
map are both 1. Two of the characteristic multipliers of any periodic solution of (|2.1|) must be 1. 
In this case all four characteristic multipliers equal 1. 

The periodic solution with the initial conditions L(0) = (p/q) 1 / 3 + O(fi), G(0) = (p/q) 1 / 3 ^ — 
e 2 ) 1 / 2 , 1(0) = n/7r, and #(0) = n g 7r depends analytically on fi at fi = and the pth return time, 
which is also the period, is given by T = 2ixp + O(fi). Therefore the entries of the Jacobian matrix 
in 1)2.2(1 will be perturbed by 0(/i) for /i > and fi sufficiently small. We determine only the 
perturbation to the entry in the upper right corner which is equal to in (|2.2|) . By (j2.1|) . 

L(T) = L(0) + /i ( T ^li(L(t),G(t),l(t),g(t))dt + 0(^ 2 ). 
Jo 
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As the solution depends analytically on \x at \i = 0, we have 

L(t) = (p/q) 1 / 3 + 0(^) 
G(t) = (p/q) 1 / 3 (l-e 2 ) 1 / 2 + 0(^ 
l{t) = run + qt/p + 0(ji) 
g(t) = rig-ir - t + O(ii) 

for < t < T. The pth return time T itself can depend upon 1(0). However, = at fi = as 
the pth return time equals 2np for all points on the L-l plane when \x = 0. Therefore 

<)U I ' " ^ ^u((p/q) l, \ (p/q) 1/ \l-e 2 ) 1 '\n l TT + qt/p, n g ir-t)dt + 0(1?). 



dl(0) 



ft 



Define C(e,p,q) such that the integral in the above equation is given by -C(e,p,q)p 1 / 3 /(6Trq 4 / 3 ). 
Then the Jacobian matrix for the pth. return map is given by 



It follows that the characteristic multipliers are 1 ± y^C(e,p, qffj, + O(fi) in the limit ix — > 0. If 
C(e,p, q) > 0, the periodic solutions in the family are hyperbolic for ti > and ji sufficiently small. 
If C(e,p,q) < 0, those periodic solutions are elliptic. 
The expression for C(e,p,q), namely 

C(e»P»?) = -^r I MM*) 1 ' 3 ' (p/5) 1/3 (l-e 2 ) 1/2 ,^ + ^/P,^-t)dt (2.3) 

can be cast into a simpler form. If Q has the same arguments as in (|2.3|) . it is periodic in t with 
period 2-rrp. Using the identity ^ = — f2; 5 , the integrand in (|2.3|) can be changed to VLi g . Use 



of the identity ^f- = ^i g — Q gg gives the expression 

C(e,p,q) = f^Qggfa/q) 1 ' 3 , (p/q) l/3 (l - e 2 ) 1 ' 2 , nm + qt/p, n g vr - t)dt. (2.4) 



7 2 / 3 .... 

The advantage of Q2.4j) becomes evident when the expression for Vt in terms of r and 9 is 

considered. From (|1.4j) and l|1.5|) . we see that r depends on L, G, and I but not on g and that 
8 = v — g, where v too depends upon L, G, and Z but not on g. Therefore 

6vrp 5 / 3 f 2np / 1 \ cos0 , 

where Ai = (1 + r 2 — 2r cos #) 1//2 . The variables r and can be obtained in terms of L,G,l,g 
using (|1.4f) and (|1.5j) . For the family of periodic solutions determined by p, q, e, rii, and n s , the 
variables L, G, Z, 5 depend upon t as indicated by the arguments of and S9 in (|2.cS|) and (|2.4|) . 
respectively. 

For the family of periodic solutions to exist, we required e £ (0, 1) to be such that Ai > 
everywhere along the unperturbed orbit at 11 = 0. The integrand of (|2.5[) can be differentiated 
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with respect to e in a complex neighborhood of the value of e for < t < 2np using formulas, since 
Ai > 0, and is therefore analytic in e. That C(e,p, q) is also analytic in e follows from a standard 
argument (see [7]). If p/q ^ 1/1, then Ai > along the unperturbed orbit for e sufficiently small. 
Therefore C(e,p, q) is analytic at e = if p/q j£ 1/1. The expression 1)2. 5|) can be used to investigate 
the expansion of C(e,p,q) in powers of e assuming p/q ^ 1/1. 

2.2 Expansion of C(e,p, q) 

Expression (|2.5j) for C(e,p,q) is an integral over an elliptic solution of the two-body problem, 
obtained by setting fi = in (|2.1|) . that satisfies Ai > at all points on the solution. By (|2.1j) and 
1)1. 5j) . / = qt/p + niiT and I = E — esinF. Change the variable of integration in 1)2.5)1 to F using 
dt = (p/g) 1 ' 3 r <iF, and then to F, which is defined by qF = E, to get 

, 6np 2 f 27Tq /r\ cosO n 2 f** / r \ cos# , n . 

C(e,p,q) = — / — + dE = -6-K P 2 — + dF. 2.6) 

q J \A 1 Jee r F J \A 1 Jee r 

Define Ci(e,p,q) and C2(e,p,q) by 

f 2jT / r \ f 2jT ens 
Ci(e,p,g)=/ ( — ) and C 2 (e,p,q) = dF. (2.7) 

Then C^e,p,q) = —Qirp 2 (Ci(e,p,q) + C^{e, p, q)) The dependence of r and on E or F is given by 
r = (p/q) 2 / 3 (1 — e cos F) and 9 = v + n s 7r — t. The variables t and f can be obtained as functions of 
F or F using I = qt/p + nitr, (|1.4|) . and (|1.5|) . The integrands of both integrals in 1)2 .7j) are periodic 
in F with period 2tt, and as Ai > 0, are analytic in F for < F < 2-7T. 

Let z = exp(iF). Define /? by e = 2/3/(1 + /3 2 ). Standard formulas that connect the true and 
mean anomalies with the eccentric anomaly (see [2|) imply the following: 

expzi/ = z q (l - (3z- q )(l - f3z q )- 1 

— (z~ q — z q ) J . (2.8) 

The variable r and exp(m#), n £ Z, too can be expressed in terms of z. 

r = {p/qf\l + (3 2 )-\l - [3z- q )(l - Pz q ) 
exp(m0) = (-l) nn « exp(i7rnn lP /q)z n( - q ~ p \l - [5z- q ) n {l - [3z q )~ n exp(^-(z q - z~ 9 )). (2.9) 

If the integrands in (|2.7|) are expressed in terms of z, they are analytic on the circle \z\ = 1. 
Therefore the integrands can be expanded in Laurent series in z. Only the constant terms in the 
Laurent series contribute to C\(e,p, q) and Cii&iP-, <?)• For later use, we record the following Laurent 
series: 

oo 

exp^(z q - z- q )) = £ Jk{enp/q)z kq . (2.10) 

k=— oo 

The Bessel functions are entire functions. We need the first term in the power series expansion 

k 7 

of Jk(x) about x = 0, which is for k > 0, but no others. For k > 0, J_fc(x) = (— 1) k Jj e {x). 
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When we consider C2(e,p,q), defined by (|2.7|) . and use a Laurent series of its integrand in z, 
we find that C2(e,p,q) = except when q = 1. Since cosO/r = (exp(i6) + exp(— i9))/(2r), we may 
use expressions for exp(in9), n = ±1, and r from (|2.9j) to get 

= f-) — ~^-^—(exp(mnip/q)z q ~' p F(z) + exp(— i7rnip/q)z p ~ q F(l/z)), 
r \pJ 2 

where F(z) = (-l) n » (1 - f3z q )- 2 exp(((ep)/(2q))(z q - z- q )). The Laurent series of F(z) only 
has powers of z q . Since p and q are relatively prime, none of these powers can equal z p ~ q or z q ~ p , 
except when q = 1. Therefore the Laurent series of cos6/r has no constant term if q ^ 1 and 
C 2 (e,p,q)=0. 

If q = 1, we can use the Laurent series (|2.10j) with n = 1 and the binomial series of (1 — [3z q )~ 2 
to obtain the Laurent series of F{z) and of cos9/r. The Laurent series of cosO/r implies 
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C 2 (e,p, q) = (-!)»» (-!)»«* 2 7 ri^(j p _ 1 (e^) + 2J p _ 2 (ep)/3 + 3J p _ 3 M/3 2 + ■ ■ ■ ) (2.11 



P 



for g = 1. Note that C 2 (e,p,q) = 0(e p ~ 1 ) as e — > 0, when q = 1. 

We will prove that the power series of Ci(e,p,q), defined by (|2.7|) . has no terms lower than 
elp~<7l We will also determine the coefficient of the e' p_(? i term. In the ensuing analysis, we assume 
p < q. The p > q case is similar. We obtain terms in the Laurent series of the integrand (3^) ee in 
two steps. For the first step, consider 



/(a) ^ 2 (vi + a 2 -2acosJ' 



This quantity 1(a) would equal the integrand if a = r. However, we assume a to be fixed in the 
range (0, 1). It will be set equal to (p/q) 2 ^ later. The Fourier expansion of 1(a) is 



^ 00 

1(a) = - -n 2 ab n (a)exp(in6), (2.12) 



where b n (a) are the Laplace coefficients. The Laplace coefficients are hyper geometric functions 
whose series converge for \a\ < 1. They and their derivatives satisfy a number of identities. See jS]. 

In the second step, we let a = r in 1(a) to make it equal to the integrand (5-J00 in (|2.7|) . Let 
f(x) be analytic in x at x = xq. Then f((l + 5)xq) has a convergent power series expansion in 5. 
This power series can be conveniently represented as (1 + 5) D f(xo), with the understanding that 
D stands for the differential operator x-^ and that (1 + 5) D is expanded according to the binomial 
formula [3]- Note that r = (p/q) 2 ^ 3 (l — ecosE). Let D be the differential operator c^5a - Since 
cosE = (z~ q + z q )/2, we get 

D a = (1 - ecosE) D = (1 + (3 2 )~ D (1 - (3z~ q ) D (l - (3z q ) D , (2.13) 

where the first equality defines the operator D a . By ()2.12j) and (|2. 13|) . the integrand in the definition 
(|2.7|) of Ci(e,p, q) can be represented as 



r 



- ^2 -n 2 D a ab n (a)exp(in6), (2.14) 
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which is to be evaluated at a = (p/q) 2 ^ 3 - We will use (j2.9j) and (|2.13|) to substitute expressions in 
terms of z for exp(in#) and D a , respectively. 

Define X n (A,B,C) = (1 + (3 2 ) A (1 - Pz~ q ) B (l - (3z q ) c exp((enp/(2q))(z q - z~ q )). The formal 
Laurent series of X n (A, B, C) is gotten by multiplying the binomial expansions of (1 — (3z~ q ) B and 
(1 — (3z q ) c with the series in ()2.1(J|) . Each of these three series is in terms of integral powers of z q 
and so is the resulting series for X n (A, B, C). In each of the three series the coefficient of z kq is a 
power series in e with the lowest power of e with possibly nonzero coefficient being e' fe ' and so it is 
for the series expansion of X n (A, B,C). 

The nth term in the summation in (|2.14l) is given by 

-n 2 D a ab n (a)exp(in9) = -(-\) nn 9 exp(iTmn lP /q)n 2 z n{q - p) X n (-D, D+n, D-n)(ab n (a)). (2.15) 

For the Laurent series of the term in (|2. 15f) above to have a nonzero constant term, n must be 
an integer multiple of q with n / 0. Let n = n^q. Then that constant term is ±n 2 times the 
coefficient of z n °( p ~ q ^ q in the Laurent series of X n (—D, D + n, D — n). By the previous paragraph, 
the lowest term with possibly nonzero coefficient in the power series expansion of that constant 
term is e' n °( p ~' ? )l As C\(e,p,q) is obtained in ()2.7j) by integrating the sum in ()2. 14|) . only the 
constant term in the Laurent series in z of each term of (|2.14|) makes a contribution to Ci(e,p, q). 
Thus the lowest possibly nonzero term in the power series of C(e,p,q) is e' p_l? ' and the only terms 
of (|2.14|) that contribute to the coefficient of e' p-<? ! are n = q and n = —q. 

Thus we have proved that the coefficients of all powers of e less than \p — q\ in the power series 
of C\(e,p, q) about e = are zero, if p < q. The coefficient of el p_9 l is equal to 

_ w (^ M)W (g(« r )_^_) (aMa)) 

evaluated at a = (p/q) 2 ^ 3 - By definition D is the operator a^, and b q is a hypergeometric function 
whose series converges for |a| < 1 as mentioned earlier. 

If p > q, a similar analysis proves that all powers of e less than \p — q\ in the power series of 
Ci(e,p, q) about e = are zero. The coefficient of e^ p ~ q ^ in that power series is equal to 




(2.17) 



evaluated at a = (q/p) 2/s . That C(e,p,q) = 0{e\ p ~ q \) for p/q ^ 1/1 is now clear from (ETKj). (t2~71l . 
(|2.11|) . and the power series of C\{e,p,q) given by (|2,l(i|) and (|2.17|) . 

The two families of periodic solutions of H2.1|) with L(0) = (p/q) 1 / 3 + O(fi) and G(0) = 
(»/(?) 1//3 (l — e 2 ) 1 / 2 + O(ijl) are given by 1(0) = niir and g(0) = n g ir, with n g = and n/ either 
or 1 if p is odd, and with n\ = and n g either or 1 if p is even. Since p and q are relatively 
prime, q must be odd if p is even. By (JOT!) . (pHfij) . and (f2~TTl) . the coefficients of e |p_91 in the 
power series of C(e,p,q) for the two families sum to zero. Of course, there is still the possibility 
that both these coefficients are zero. This possibility can be eliminated by evaluating (|2.16j) if p < q 
and (|2~T7|) if p > q, and also (|2~TT|) if q = 1. 

2.3 Statement of results 

The theorem below summarizes the results of this section. 
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Figure 3: The plots use p/q = 1/3 and p/q = 2/7, respectively 



Theorem 2.1. There exists a p-dependent family of periodic solutions of Hamilton's equations 
with the initial conditions L(0) = (p/q) l/3 + 0(p), G(0) = (p/q) 1/3 (l - e 2 ) 1 / 2 , 1(0) = mn, 
and g(0) = n g Tr, where p and q are relatively prime positive integers, < e < 1, and ni and n g 
are either or 1, provided the solution at p = does not collide with the orbit of the primary of 
mass p. The existence holds for p sufficiently small, the dependence on p is analytic, and there is 
only one such family as proved in 0J/ ; J3J/. Given p, q, and e there are four possible choices for ni 
and n g , but only two of these give rise to distinct families. Two of the characteristic multipliers 
are equal to 1 ± y / C(e,p, q)p + O(p) in the limit p — ► 0, where C(e,p,q) is given by (|2.5[) . The 
quantity C(e,p,q) is analytic at e = and C(e,p,q) = 0(e' p " 9 ') as e — > 0, if p/q ^ 1/1. The 
coefficient in front o/e' p_l3 ' in the power series of C(e,p,q) is given by (|2.11|) . (|2.16|) . and (|2.17() . 
The coefficients for the two distinct families obtained for given p, q, and e sum to zero. 

The numerical evaluation of C(e,p, q) defined by 1)2. 5 Jl can be carried out with great accuracy at 
little expense. Figure El graphs C(e,p, q) against e with two choices of p/q. In the graphs, C(e,p, q) 
has opposite signs for the two distinct families. For p > and p small, if the periodic solutions 
of one family are hyperbolic, the periodic solutions of the other family are elliptic in either plot of 
Figure |31 This corresponds to the situation depicted in Figured 

3 Retrograde periodic motions 

The convention for representing retrograde motions in the restricted three body problem depends 
upon the choice of variables. We will represent retrograde motions using negative L and G. The 
variation of L,G,l,g with t is again obtained from Hamilton's equations (|2.1|h The dependence 
of r and 8 on L,G,l,g is given by (|1.4|) and (|1.5|) without a change. For the choice of variables 
to be valid, the condition < e < 1 must hold. If the angles E, v, and g are interpreted in the 
inertial frame, say with p = 0, they must all be measured in the counterclockwise sense. However, 
for retrograde motion, unlike for direct motion, the angles E and v decrease with t. 

There exists a //-dependent family of retrograde periodic solutions of (|2,lj) with initial conditions 
L(0) = -(p/q) 113 + 0(p), G(0) = ~(p/q) 1/3 (l - e 2 ) 1 / 2 , 1(0) = n,vr, and g(0) = n g 7r, where p and q 
are relatively prime positive integers, < e < 1, and n\ and n g are either or 1. The dependence 
on p is analytic and the existence holds for p > and p sufficiently small provided the solution at 
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/i = satisfies Ai > 0. Given p, q, and e, there are four possible choices for m and n g . Of these 
only two give rise to distinct families. If p is odd, the two distinct families can be obtained with 
n g = and n\ either or 1, and if p is even with n\ = and n g either or 1. 

Up to terms of order /x, the dependence of periodic solutions of the family determined by p, q, 
e, ni, and n g is given by 

L(i) = -(pA?) 1/3 + O0u) 
G{t) = -{p/q) 1 ^{l-e 2 ) 1 ' 2 + 0^) 
l(t)=rn7r-qt/p + 0(jj) 
g(t) = n g TT - t + 0(n) 

with the period being T = 2irp + O(p). To find the characteristic multipliers of these periodic 
solutions, we choose the Poincare section g = n g ir and represent it using the L-l plane. As in 
Section 2, we may deduce that 

d(L(T),l(T)) _ ( 1 0\ U - C (e,p,q)£ /3 X 



<9(L(0),Z(0)) 
The quantity C(e,p, q) is given by 




C(e,p,q)=-^^- % g {-{p/q) ll \-{p/q) l, \^ ~ e 2 ) 1/2 ,n l7 r - qt/p,n g 7T - t) dt. (3.1) 



o 



This expression for C(e,p, q) is similar in form to (|2.4|) . but the arguments of Sl gg in 1)2.41) and (|3,1|) 
are different. Use the change of variables from t to F given by I = —qt/p + niTr, I = E — esinE, 
dt = —(p/qY^rdE, and qF = E, to get 



-" ' r \ cos i 



C(e,p,q) = -6np 2 J Q ^-) $9 + — dFt (3.2) 
where Ai = (1 + r 2 — 2r cos 

^) x /2. This expression for C(e,p,q) looks identical to (|2.6[) . but the 
dependence of 6* on t, and hence on F, is quite different. The dependence of L, G, I, g on t here is as 
indicated by the arguments of tt gg in (J3.1|) and not as in ()2.4|) . The variables r and are obtained 
from L,G,l,g using (|1.4|) and ()1.5() as before. We now consider Ci(e,p,q) and C2(e,p,q) which 
are defined as in ()2.7|) but with the dependence of r and 9 on F that is indicated here. We have 
C(e,p,q) = -Q-Kp 2 {Ci(e,p,q) +C 2 {e,p,q)). 

That C(e,p,q) is analytic in e follows from the assumption Ai > as in Section 2. Here too 
C(e,p,q) is analytic at e = \ipjq / 1/1. But the power series of C(e,p,q) about e = looks 
quite different, as will be shown now. 

Let z = exp(iF). Define (3 by e = 2/3/(1 + (3 2 ) as before. We have 

exp(ii/) = z 9 (l - /3z" 9 )(l - Z?^)" 1 

— (z -9 — z q )J , (3-3) 
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which has —t in place of the t in (|2.8[) . From r = (»/(/) 1//3 (l — ecosE), 9 = v + n g 7r — i, and (|3.3j) . 
we have 

r = (p/g) 2/3 (l + (3 2 )-\l - Pz~ q )(l - (3z q ) 
exp(mfl) = eM-wnn lP /q)(-l) nn v z n(p+q \l - [5z~ q ) n {l - pz q )~ n exp(^-{ Z - q - z q )). (3.4) 

The expression for exp(in#) in (|3.4|) has a factor while the expression in ()2.9|) has the factor 

z n (i-p) , A consequence of the difference in the two expressions is that C(e,p,q) for retrograde 
motion is 0(e p+q ), and not 0(e' p_l3 '), in the limit e —* for p/q 7^ 1/1, as will be shown now. 

The analysis of C2(e,p,q) is similar to that in Section 2 with the difference that (|3,4|) and not 
(|2.9|) must be used to express cos(6?) in terms of z. We can conclude that C2(e,p,q) = if q 7^ 1 
and that C2(e,p, q) = 0(e p+l ) when q = 1. Besides, the sum of the coefficients of e p+1 in the power 
series of C(e,p, q) about e = for the two distinct families, obtained for given p, q, and e, is always 
zero. 

We consider C\(e,p, q) with the assumption p < q. The p > q case is similar and leads to the 
same conclusions. The integrand {-%-)g e can again be written as a sum of terms as in (|2.14|) and 
the operator D a is again given by (|2.13() . But the expression for the nth term in (|2.14|) is now 
different. Instead of ()2.15|) . we now have 

-n 2 D a ab n (a) exp(in#) = -(-\) nn v exp(-i7mn l p/q)z n{p+q) X n (-D, D + n,D — n)(ab n (a)), (3.5) 

where X n (A, B, C) = (1 + (3 2 ) A {1- f3 Z - q ) B {l- j3z q ) c exp(^(^- (? - z q )). The Laurent series of the 
nth term of l|2.14|) . which is displayed in (|3.5j) . has a possibly nonzero constant term only if n = noq, 
where no is an integer, and no 7^ 0. That constant term is ±n 2 times the coefficient of z n °( p+q ^ q in 
X n (—D, D + n, D — n), and therefore, if that term is expanded about e = the lowest term with 
a possibly nonzero coefficient is e' n °( p+<? )'. Thus the lowest term with a possibly nonzero coefficient 
in the power series of C\(e,p,q) about e = is e v+q and the only terms ([2.14ft which contribute to 
its coefficient are obtained by setting n = ±q in (|3.5j) . We conclude that C\(e,p,q) = 0(e p+q ) as 
e — ► for p/q ^ 1/1. 

For given p, q, and e, the two distinct families are given by n g = and n\ either or 1 if p is 
odd, and by ni = and n g either or 1 if p is even. Inspection of (|3.5|) with n = ±q leads to the 
conclusion that the sum of the coefficients of e p+q in the power series of C\(e,p, q) about e = for 
the two families is zero. 

The theorem below is about retrograde periodic motions. 

Theorem 3.1. There exists a fi- dependent family of periodic solutions of Hamilton's equations 
(t2~Tjl with the initial conditions L(0) = ~(p/q) 1/3 + O(fi), G(0) = ~(p/q) l/3 (l - e 2 ) 1/2 , 1(0) = run, 
and g(0) = n g ir, where p and q are relatively prime positive integers, < e < 1, and n[ and n g 
are either and 1, provided the solution at fi = does not collide with the orbit of the primary of 
mass 11. The existence holds for /i sufficiently small, the dependence on 11 is analytic, and there is 
only one such family. Given p, q, and e, the four possible choices for ni and n g give rise to only 
two distinct families. Two of the characteristic multipliers are equal to 1± ^JC(e,p, q)n + O(fi) in 
the limit 11 — > 0, where C(e,p,q) is given by (|3,2|) . The quantity C(e,p,q) is analytic at e = and 
C(e,p,q) = 0(e p+q ) as e — > 0, if p/q 7^ 1/1. The sum of the coefficients of e p+q in the power series 
of C(e,p, q) about e = for the two families is zero. 
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4 Periodic motions near collision 



In this section, we consider periodic motions near collision with the primary of mass Consider 
the Hamiltonian H defined by The generating function S = (— p, + £ 2 — v 2 )p x + (2£v)p y can 

be used to effect the Levi-Civita transformation from the variables p x ,p y ,x,y to the variables 
P£,Pv,£,v. The Hamiltonian becomes 



Suppose we are interested only in solutions of Hamilton's equations with H = C. Consider K 
{i 2 + v 2 ){H-C) or 



2 ((^2_^2_ 1 )2 +4 ^2) 1 /2- 

The solutions of Hamilton's equations formed using K with K = correspond to solutions of 
Hamilton's equations of H with H = C, but with time rescaled from t to r such that dt = 
(£ 2 + u 2 )dr. Hamilton's equations of K admit solutions with K = that pass through or close to 
£ = f = 0. These can be interpreted as analytic continuations of the solutions of the equations 
formed using H through the singularity of H at the location of the primary of mass 1 — [i. When 
[i = 

r>2 i 2 ^2 I 2 

* = 8^ + 2^ " ^ " 2C) " L 
The quantity G = £p„ — i^pg is a first integral of the Hamilton's equations of K shown in Q4.1|) . 
If a solution of those equations passes through £ = z/ = 0, then G = 0. In terms of p x ,p y ,x,y, 
G = 2(xp y — yp x )- Therefore the quantity denoted by G in this section is twice the angular 
momentum. 

The Hamilton-Jacobi equation of K shown in Q4.1|) was solved in jS] in the region of phase space 
with G = or G ~ 0. There is a gap in that solution when G = 0. In this section, we discuss that 
gap and make minor corrections to the solution given in [5] for G ^ 0. We follow the geometric 
approach to the construction of action-angle variables described in [2] . All references to [2] in this 
section are to the last chapter of that book. 

The generating function S = p^(r cos 6) + p v {r sin#) can be used to effect a change from the 
variables £>£,Pi/,£,f to the variables R,G,r,8. The change to the polar variables r, and the 
corresponding generalized momenta R, G is well defined only if r ^ 0, or equivalently, £ 2 + v 2 ^ 0. 
The Hamiltonian K becomes 

i (~~<2 2 2 

K= ±[R 2 + — ) +?_(_£ -2C)-l + ^(i?rsin(20) + Gcos(20)) + //- '"" 



r 



2 V ; 2 K K ' v " ^ (l + r 4 -2r 2 cos(20)) 1 /2' 



i . (72. 2 

tf = -( J R 2 + _) +_ (_G_2C)-1. (4.2) 



When [i = 



r ) 2 

We will find the action-angle variables of K at [i = using (|4,2j) . 
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Figure 4: The cycles used to obtain action-angle variable are shown above. 



The action-angle variables will be found in the region of R, G, r, 6 space or of P£,Pv,£,, v space 
where the following conditions are satisfied: 

G + 2C<0, K + 1>0, (K + l) 2 + G2 ( G + 2C ) > 0; (4.3) 

with C < 0. The necessity of these conditions will become clear shortly. We are interested in 
solutions that are near a collision with the primary of mass 1 — fx. Therefore G ~ in the region 
of interest. Further, the solutions of interest are bounded and satisfy K = 0. Therefore C < and 
K + 1 ~ 1 in the region of interest. It is obvious that all points in the region of interest meet the 
conditions 1)4.311 . 

As a first step to the derivation of action-angle variables, we visualize the set of points in 
R,G,r,9 for fixed values of K and G satisfying (|4.3|) and G / 0. Using 1)4.2)) . we deduce that 
Train < r < r max where 

2 2 2{K + l) 2 2 G 2 , AA . 

r ■ ■ + r = — — r ■ r = (4 4) 

' mm 1 ' max / q 2C) ' mm Q 2C) ' ' 

The last of the conditions ()4.3|) ensures < r min < 

r max- The condition G ^ implies < T m i n . 
For later use, we introduce variables a and e defined by 

r min = e), r 2 max =o(l + e). (4.5) 

For r within the allowed range, any value of 8 is allowed in the set. Thus the projection of the 
set looks like the annulus in Figure 0J The variable G is of course fixed. For a given point in the 
annulus in the r-0 plane, R can have two values of equal magnitude but opposite sign. On the 
boundary of the annulus, R = 0. Therefore the set of points is a torus. 

The set of points with G = and a fixed value of K, such that the conditions ()4.3)1 are met, 
includes points with £ = v = 0. The change to polar variables is invalid at those points and the 
R,G ,r,9 space cannot be used to visualize that set. By 1)4. the set of points in the P(,Pv,£,,v 
space is given by the conditions 

£ Pl/ - m = 0, ? + l ?<(K + l)/(-C), pl+pl = 8(K + l + C(e + v 2 )). 
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This set of points is also homeomorphic to a torus. It is easy to check that the differential 1-forms 
dG and dK are linearly independent on this torus. 

We return to Figure |I] and the assumption G / and find the action- angle variables. As 
explained in |2j, the action variables are functions of K and G obtained as f Rdr + Gd9, with 
the integral taken over two linearly independent cycles on the torus. Our choice of cycles is shown 
in Figure 0] The cycle with fixed r and 9 varying from to 2ir gives G as an action variable. The 
other cycle increases from r m i n to r max with R > and then decreases to r m i n with R < 0. The 
variable 9 is fixed on this cycle. The action variable L* given by this cycle is 



1 fVmax , (~<2 \ 1/2 

8^ + 8- — + 4r 2 (G + 2C)) dr. 



7r L. . \ r 



Substitute u = r 2 to get 



L* = k-C - 2C)V 2 f ^ (WH.)-)(»-.(l- e ))) 1/ ' fc 



Ja(l-e) 

Evaluate the integral and use (|4.4j) and (|4.5|) to get 



= (-G-2C)V* " T" (46) 



Instead of L*, we use L defined by 



iT + 1 

L = (-G-2C)^ (47) 



as the other action variable. 

The angle variables parametrize the surface of the torus obtained in R, G, r, 9 space by fixing 
the action variables L and G. As proved in [2j, the angle variables can be derived from the following 
generating function: 



S(r,9,L,G) = / Gd9 + Rdr. 

" fmAn < 



,0 

In this expression, it! must be expressed in terms of L, G, and r using (|4.2|) and (|4.7|) . From this 
expression, it might seem that S(r, 9, L, G) is a function on the annulus shown in Figure 0] for fixed 
L and G. It is actually a function on the torus as R can be either positive or negative. 



S(r,9,L,G) = G9± f I 8L(-G - 2C) 1/2 - ^ + 4r 2 (G + 2C)\ dr, 

^ 'min y / 

where the sign is + if R > and — otherwise. The generating function S(r, 9, L, G) is a multiple 
valued function on the surface of the torus as there are many non-homotopic paths from the base 
point R = 0, r = r m j n , 9 = to any point on the torus. In the calculations below, we consider the 
paths for which both r and 9 increase monotonically. 

The value of one of the angle variables on a point on the torus is equal to the value of M at 
that point on the torus as proved in 



^- = ±4(-G - 2C) 1/2 / 8L(-G - 2C) 1/2 - ^ + 4r 2 (G + 2C) 
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where the sign is + if R > at the point on the torus and — otherwise. Change variable to u = r 2 
to get 

dS f r2 du 



dL Ja(l-e) (( (1 +e ) - U )(u - (1 - e ))) 1/2 ' 

Change variable to I defined by u = a(l — ecosZ) with < Z < ir to get ^ = ±Z. If we adopt the 
convention 0</<7rifi?>0 and —tt < I < if R < 0, we get 

where Z is an angular variable that is measured modulo 27r. 

The value of the other angle variable at a point on the torus is equal to the value of at that 
point on the torus. 



^ = 0±±jT L 4 L(-G-2C7)-^_^ + 4r 2j/ 8L( _ G _ 2C) i/2_^ + 4r 2 (G + 2C) j ' dr 



, -1/2 
2 \ 



where the sign is + if R > at the point on the torus and — otherwise. Change variable to u = r 2 
to get 

8S „ 1 T 2 (-2L(-G-2C)" 1 / 2 -G/u + 2n) , 

/- 1 ' ■ du. 



dG 4(-G-2C)V2y a(1 _ e) ((a(1 + e) _ n)(n _ a(1 _ e)) i/2 

Change variable to Z defined by u = a(l — ecosZ) with the same convention as in (|4.8|) . By (|4.4|) . 
dUSD, and @2J, 

-2L(-G - 2C) 1/2 - G/u + 2u = 7 - 2ae cos Z. 

a(l — ecosZ) 



Denote §^ by g to get 



G /•' dZ (L 2 — G 2 /4) 1 / 2 , 

' smZ, (4.c 



4L 7 (1-ecosZ) 2(-G-2C) 
where e = (l — Ifr) 1 ^ 2 - 

If L*, defined by l|4.6jl . and G are used as the action variables, the corresponding angle variables 
are Z and g±l/2 — the sign is + if G > and — if G < 0. It is easily verified that along each of the 
two cycles depicted in Figure |1] one of these angle variables increases by 2ir while the other doesn't 
change. This must be the case as proved in [2] and the verification is a check of the correctness of 
our derivation of action-angle variables. 

The quantities denoted by G and g in jS] are twice and half of the quantities denoted by G 
and g in this section, respectively. The factor 2 that appears in the coefficient of sinZ in (19) of [0] 
must be moved to the denominator. Further, in (20) of 6 (1 — e 2 ) 1 / 2 must be replaced by G/L. 
Although | G/L | = (1 - e 2 ) 1/2 in the notation of the sign of G is significant. 

When G = the change to polar variables is not valid as we have shown in this section. 
Therefore the derivation of action-angle variables using a change to polar variables given in [2j is 
also not valid when G = 0. 
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The existence of /i-dependent families of periodic solutions of Hamilton's equations of K with 
the initial conditions L(0) = {p/q) 1/3 + 0(jj), G(0) = -2C - (q/p) 2/3 + 0(p), 1(0) = 0, and 
g(0) = was proved in ^Uj. The solutions in these families depend analytically on fx and K = 
at all points along the solutions. The methods of Sections 2 and 3 can be used to carry out linear 
stability analyses of these periodic solutions. However, a new derivation of action-angle variables 
is necessary for the crucial G = case. 
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